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This paper presents a calculation of the neutron cross-sections in solid materials (used in practical 
neutron sources) with a large coherent scattering contribution. In particular, the dynamic struc- 
ture function S(Q, to) of polycrystalline ortho-D2 is evaluated using a Monte-Carlo calculation that 
performs an average over scattering angles relative to crystal axes in random orientations. This 
method uses an analytical dispersion function with force constants derived from neutron scattering 
data of single crystal D2 in the framework of an axially symmetric force tensor. The resulting two 
dimensional map of S(Q, uj) captures details of the phonon branches as well as the molecular ro- 
tations, that can be compared directly to data from inelastic neutron scattering on polycrystalline 
D2. This high resolution information is used to calculate the absolute cross-sections of production 
and upscattering loss of ultracold neutron (UCN). The resulting scattering cross-sections are signif- 
icantly different, especially for UCN upscattering, from the previous predictions using the approach 
centered on the incoherent approximation. 
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I. INTRODUCTION 



Neutron moderation is a well established technique, 
developed for applications in fission reactors [l| and nu- 
clear weapons. Interest in this topic has been revived 
with the development of new neutron scattering facilities 
for investigation of materials, new sources of ultracold 
neutrons for fundamental physics, as well as new deep- 
underground laboratories which require extensive shield- 
ing. In the latter case, neutrons generated in the sur- 
rounding rock via the spallation process initiated by high 
energy cosmic rays are a major source of background. 
On the pure academic front, neutron moderation can be 
put in a broader context of phase space compression, de- 
veloped for several different types of particle beams in 
various areas of physics, from stochastic electron cooling 
in high energy particle accelerators on Tera-eV scales to 
laser cooling of neutral atoms to form Bose-Einstein con- 
densates on pico-eV scales. The phase space compression 
of neutrons is especially challenging. Because the electro- 
magnetic interaction of the neutron is relatively weak due 
to its zero charge and small magnetic moment, the usual 
mechanisms of phase space compression for charged par- 
ticles do not apply. Instead, the nuclear force couples 
neutrons to individual nuclei in the interacting medium. 
Upon scattering, the nuclear recoil leads to rapid mod- 
eration for fast neutrons, but stops to be effective for 
neutrons with energy below a few tens of milli-eV. For 
neutrons of particular interests to fundamental research, 
such as cold, very-cold, and ultra-cold neutrons (with 
energies ranging from milli-eV to nano-eV), the only effi- 
cient way to date to increase neutron phase space density 
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is through inelastic collisions in a cold medium, so-called 
neutron moderator. 

To study the process of neutron cooling and estimate 
the neutron fiux for the broad range of applications men- 
tioned above, standard simulation codes are available in- 
cluding the widely-used Monte Carlo N-Particlc Trans- 
port Code (MCNP)0. The code was developed in the 
1960s at Los Alamos National Laboratory, and has since 
been extensively bench-marked with experimental data; 
the physics models and scattering kernels used have been 
refined with every new release. In MCNP, moderations 
of neutrons in a medium are modeled by the Boltzmann 
equation that details the evolution of the phase space dis- 
tribution function of neutrons. The Boltzmann equation 
includes neutron diffusion, neutron absorption loss, and 
inelastic processes in which neutrons downscatter with 
energy loss or upscatter with energy boost. In this for- 
malism, the scattering kernel determines the probabil- 
ity for inelastic processes. For fast neutrons (with en- 
ergies from cV to MeV), nuclear resonances strongly in- 
fiucncc the scattering and absorption processes. Away 
from the resonances, the neutron energy loss through nu- 
clear recoil is calculated with a simple free gas model, in 
which the time scale of collisions is so short that nuclei 
are treated as free. For neutrons with wavelength larger 
than the separation between scattering nuclei, however, 
thermal neutron scattering kernels that include coher- 
ent diffraction and collective excitations in bulk materi- 
als have only been incorporated for limited materials [3| . 
Whereas nuclear cross-sections for fast neutron reactions 
are detailed for almost every stable isotope, there ex- 
ist only a handful of materials with extensive treatments 
that include condensed matter effects beyond the simple 
nuclear recoil in the free gas model. 

Hydrogen and deuterium in liquid and solid forms are 
of particular interest, as they are the materials of choice 
to construct cold neutron moderators used in research 



reactors and spallation neutron sources. The large nu- 
clear cross-section, together with the maximum energy 
loss in recoil kinematics of n-^H scattering, leads to rapid 
cool-down of incident neutrons. However, the -'^H(n,7)^H 
reaction results in significant loss of the population of 
the already moderated, low-energy neutrons. Neutron 
absorption loss in ^H (i.e., D) is much less of an issue, so 
even though the cross-section is an order of magnitude 
smaller, D compares favorably to H as a good neutron 
moderator. Since the ground state of hydrogen forms a 
diatomic molecule, the scattered wavefunctions of the in- 
cident neutron off the two identical nuclei within a single 
molecule interfere coherently. Together with the spin- 
dependent nuclear force, the inter-molecular interference 
leads to different scattering amplitudes between states 
with distinct nuclear spin. Furthermore, spin statistics 
modify the molecular form factor, giving rise to selec- 
tion rules governing the coupling to discrete sets of rota- 
tional states [J|. In fact, applying Fermi-Dirac statistics 
in the diatomic H2 molecule was the key to solve the puz- 
zle of discrepancies between measured values of nuclear 
cross-section of neutron-proton scattering and theoreti- 
cal predictions using the nuclear bound state of D [a, , 
during the early construction of nuclear structure theory. 
To this end, most cross-section evaluations have properly 
included this coherent spin effect 0] . 

On the other hand, all reported works (including 
NJOY 01, d], [1], M) that evaluate inelastic cross- 
sections in the thermal and cold regimes for use in MCNP 
employ the incoherent approximation (lA) that assumes 
negligible interference between wavefunctions scattered 
from different lattice sites in the intra-molecular scale. 
In the lA, the amplitude of coherent scattering is eval- 
uated using the same algorithm used in evaluating the 
incoherent scattering. The use of lA is popular because 
the inelastic incoherent cross-sections can be readily cal- 
culated via standard phonon expansions, only requiring 
limited information on density of states (or the frequency 
distribution, as referred in some literature) of the col- 
lective excitations as an input. While lA works quite 
well for moderators consisting of hydrogenous materi- 
als (erf '^=1. 7583 b, o-f ^==80.27 b), in the case of D2, 
the contribution of the coherent scattering is not small 
(CTfj°''=5.592 b, cr™'==2.050 b) and thus the adoption of 
lA should be examined carefully. Especially for low en- 
ergy neutrons with wavelengths comparable to and larger 
than the lattice constants, the use of lA seems particu- 
larly deficient. To address the shortcomings of this as- 
sumption and all previous calculations that invoked this 
assumption, we have developed the first full treatment of 
neutron scattering for neutrons with wavelength larger 
than angstroms. This calculation includes both the co- 
herent and the incoherent scattering which excite single 
and multiple phonons, as well as the processes which ex- 
cite rotational transitions. The elastic Bragg scattering 
is also included. In this paper, we will show that this 
full model, when applied to solid D2, significantly modi- 
fies the cross-sections of ultracold neutrons (UCN), which 



arc free neutrons with E/fc^ < 3 niK. where ks is the 
Boltzmann constant. 

The scattering process responsible for UCN production 
differs from the thermalization process in the thermal 
and cold neutron moderators, in which multiple scatter- 
ings take place while the energy spectrum evolves con- 
tinuously until the neutron gas establishes thermal equi- 
librium. In contrast, superthermal UCN production oc- 
curs as a single scattering event during which the inci- 
dent neutron comes to a near-full stop, giving up its en- 
ergy and momentum to a matching quasi-particle created 
in the interacting media. Even though the phase space 
of this inelastic scattering process is limited, the cross- 
section is non-negligible, provided that the interacting 
medium has excitations that coincide with the energy 
and momentum of incident cold neutrons. Subsequent 
scattering of the down-converted UCN (mostly upscat- 
terings) is suppressed inside a superthermal UCN source 
by simply reducing the thermal population of the quasi- 
particles. Cooling the moderator effectively pumps out 
these quasi-particles. For the typical application of neu- 
tron thermalization down to a few milli-eV, scattering 
kernels evaluated based on lA give reasonable predic- 
tions of the neutron yields. However, many criteria of 
lA break down when applied to the extreme low-energy 
regime (E < 350 neV) of UCN physics. Rc-evaluating 
the UCN cross-sections in solid ortho-D2 using the full 
model shows that UCN cross-sections differ significantly 
from previous estimates using the lA approach [8|. 



II. DYNAMICS IN SOLID D2 

Calculation of the coherent inelastic scattering in solid 
D2 requires detailed knowledge of the dynamics and the 
dispersion of the collective energy excitations. Fortu- 
nately, the dynamics of D2 in a molecular lattice were 
studied extensively in the early days of neutron scattering 
experiments |lll-[l9|. Due to the large lattice constants 
resulting from the large zero-point motion, molecular ro- 
tation remains free in the solid matrix and the trans- 
lational degrees of freedom can be decoupled and the 
coupling strength estimated using an isotropic Lennard- 
Jones potential. The crystal is simple enough that the 
dispersion energy and the polarization vector of distinct 
phonon branches can be calculated with a tensor-force 
model, without the use of molecular dynamic simula- 
tions that require intensive numerical computations. De- 
wames et al. [20l - |23 | developed a Born-von Karman 
model, using an axially symmetric (A-S) two-body po- 
tential including up to the third nearest-neighbors, for 
systems with the hexagonal close packed (HCP) struc- 
ture. This dynamical matrix model was applied success- 
fully to HCP terbiu m |23ll , and recently extended to HCP 
Tb, Sc, Ti and Co [24l . Using the force constants mea- 
sured by Nielsen |18| in the A-S model framework, we re- 
constructed the full energy dispersion that is a function 
of the three dimensional momentum vector of phonons 



propagating in solid ortho-D2. This angular-dependent 
dispersion function is then used in a Monte-Carlo code to 
construct the inelastic coherent scattering cross-section 
for polycrystalline solid D2. For each scattering event 
sampled by the Monte-Carlo, the scattering kinematics is 
determined and the coherent scattering amplitude is cal- 
culated following the standard textbook formalism [25| . 
The major difference in evaluating the coherent and in- 
coherent scattering cross-sections lies in the kinematic 
conditions of momentum conservation. For the incoher- 
ent process, the scattered wavefunctions are not summed 
coherently as the phase coherence is destroyed by fluc- 
tuations of the scattering sites. This leads to a relax- 
ation of the momentum conservation law typically ap- 
plied to two-particle systems. Instead the approach of 
lA uses the density of states of the phonon modes to 
estimate the relative contribution of the scattering am- 
plitude for phonons of different energy. The energy of 
phonon equates to the energy transfer of the scattered 
neutron a result of energy conservation. As a conse- 
quence, cross-sections evaluated using lA has a diffuse 
smooth Q dependence without any localized enhanced 
intensity, as would have been expected in coherent scat- 
terings. 

Solid hydrogen and deuterium are considered quan- 
tum solids due to the light mass of individual molecules. 
The quantum nature is determined by its thermal wave- 
length, defined as At = 2tt M/^/iWT [HI. For D2, it 
is 4.6 A at 5 K and 3 A at 18 K. For a light molecule 
(and atom), its large thermal wavelength exdending be- 
yond the inter-molecular spacing causes direct overlap 
of the wavefunctions of adjacent molecules. This could 
lead to quantum inteference. Provided that the adja- 
cent molecules are identical particles, the Bose-Einstein 
statistics could significantly modify the collective behav- 
ior of the molecules, as in the case of superfluid helium. 
Here, the ortho-state of D2 molecule has possible nuclear 
spin of and 2 and distinct sub-states of |/, m) = (0,0), 
(2,2), (2,1), (2,0), (2,-1), (2,-2) that are equally popu- 
lated. Given that the seperation between molecules on 
the HCP lattice is d = 3.607A, the mean seperation of 
particles with identical states is ^y/& x d = 6.55 A, which 
is larger than the thermal wavelength by ^ 40%. This 
suggests that at temperatures higher than 5 K, solid D2 
can be treated as a classical system to the first order ap- 
proximation. Many ongoing research [ij, |26| are looking 
for the evidence of quantum interference in H2 , HD and 
possibly D2. 

In a classical system, the double differential cross- 
section of neutron scattering is conventionally expressed 
as: 



to the scattering kernel used in the neutron-moderation 
Boltzmann equation, in which the rate of neutron mod- 
eration is directly proportional to the dynamic structure 
function of the interacting medium. For the remaining of 
the paper, we examine neutron scattering in a classical 
solid lattice, where the scattering can be described as 
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The scattering length b is determined by the strength 
of interaction between the neutron and the scattering 
site. The Debye- Waller factor M^ is a result of the de- 
localization of the atoms on the lattice sites due to ther- 
mal motion (as well as zero-point motion), smearing out 
the sharpness of constructive interferences of the scat- 
tered wavefunctions. On a lattice with simple cubic sym- 
metry, at low temperatures, W reduces to 
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where M is the mass of the lattice site. For solid D2, 
M is the total mass of the D2 molecule. The Debye 
temperature ujd for solid D2 is around 110 K. 
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where k/ko arises from the ratio of phase space before 
and after scattering, and the dynamic structure function, 
S{Q,uj), contains all the detailed physics of the neutron 
interactions. The dynamic structure function is related 



FIG. 1: The dynamic structure function of solid ortho-D2 at 5 
K, calculated using incoherent approximation (lA), presented 
as a contour plot. The color bar on the left indicates the 
intensity of the structure function. Both 1-phonon and 2- 
phonon processes are included. The solid curve is the free 
neutron energy dispersion curve, which outlines the kinematic 
condition for UCN production. 



The first term in Eq. ^ is the coherent elastic scatter- 
ing, during which the energy transfer of the neutron is 
zero but the momentum transfer can be absorbed by the 
whole lattice if the Bragg condition {Q = 2ksm{6/2) ~ 
G) is satisfied. The next term describes the inelastic 
scattering where the neutron energy is released to cre- 
ate a single phonon in the j*'' mode with energy u!j{q) 
and polarization vector Cj. The occupation number of 
the created phonon is nj{q), which satisfies the Bose- 
Einstein statistics. The delta functions enforce energy 
and momentum conservation. Any change in momentum 
of the incident neutron during collision is picked up by 
a phonon of momentum q. The total momentum con- 
servation is adjusted by the appropriate inverse lattice 
vector G if the change of momentum is beyond the first 
Brillouin zone. The third term leads to energy and mo- 
mentum gain by the neutron through absorption of a sin- 
gle phonon. Higher-order terms describe multi-phonon 
processes, with rapidly decreasing probabilities for small 
momentum transfer Q. Finally, with D2 molecules ar- 
ranged in a solid lattice with lattice vector 0,1,0,2,0,^, the 
inverse lattice vectors G along the principle axes are 
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where V = ai . [0,2 x as) is the volume of the unit cell. 



III. ELASTIC SCATTERING AND FORM 
FACTORS 

Even though the clastic scattering is well understood, 
comparing the experimental data of cross-sections with 
the calculation allows for consistency checks on the 
occasionally non-trivial molecular form factors. This 
is particularly useful when calculating cross-sections in 
solid oxygen, where the spin form factor is less well- 
known than the form factors of atomic nuclei, with ad- 
ditional complications due to preferential alignments of 
the molecular axis relative to the crystal axes. For D2, 
the situation is somewhat simpler because the rotational 
motion remains free in the solid lattice, and the molecu- 
lar wavefunction can be described by spherical harmonics 
without preferred orientations [Jl . The total cross-section 
for elastic scattering (integrating the first term in Eq.®) 
can be separated into coherent Bragg part and incoherent 
diffuse part [2l,[23|: 
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where d/i^i is the inter-planar spacing between the [hkl] 
planes in the lattice space, i.e.. 
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The form factor F within the unit cell differs depend- 
ing on whether the coherence of the scattered wave is 
preserved after scattering. For coherent processes, 
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In a non-Bravais crystal, the number of atoms in the unit 
cell, Nfjasis, is larger than 1. 

On each lattice site in the D2 crystal resides a diatomic 
molecule, the additional molecular form factor of which 
significantly modifies the scattering length. For D2, the 
bosonic nature of the two identical deuteron nuclei leads 
to unique couplings to a set of molecular rotational states 
of either even or odd symmetry, so that the symmetry 
of the whole system is preserved under permutation of 
identical particles. It follows that there exist two differ- 
ent types of D2 molecules, i.e., ortho-D2 with even nu- 
clear spin (/=0, 2) coupling to symmetric molecular ro- 
tational wavcfunctions (with orbital angular momentum 
J=0,2,4,...), and para-D2 with odd nuclear spin (/=1) 
coupling to anti-symmetric molecular wavefunctions (of 
J=l,3,5,..). By convention the ortho-state designates the 
group with the highest spin multiplicity. On account of 
the spin-dependence of the nuclear force, neutron scat- 
tering lengths of ortho and para-D2 are different. In 
addition, when interacting with one species of diatomic 
molecule, say ortho-D2 molecules (with 1—0,2 and J=0 
in the ground state), the neutron scattering length also 
depends on the final state of the target molecule [J], i.e.. 
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where acoh ~ 6.67 fm and Uinc 
ular form factors are: 
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= 4.04 fm, and the molec- 

A02 = 2j2 



(12) 
with jl as the spherical Besscl function of order I. The 
separation of the two deuteron nuclei is a = 0.71 A. All 
Clebsch-Gordan coefficients used in the above equations 
arc 1. The J = — >■ 1 process resulting in the ortho- 
to para- transition requires that the neutron to transfer 
energy equal to the rotational transition energy, and thus 
it is only present in the inelastic scattering. 

The energy dependence of the total elastic scattering 
cross-section of polycrystalline samples is determined by 
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around the direction of the c-axis. The reduced scattering 
cross-section could also be explained by multiple scatter- 
ing processes, in which some neutrons that are already 
scattered once scatter again back into the beam and thus 
enhances the neutron transmission. The presence of mul- 
tiple scattering leads to a reduction in the cross-section 
if the analysis algorithm does not include a correction for 
multiple scattering. Due to this complication, the data of 
neutron transmission, which was collected along with the 
data of UCN production, cannot be used to distinguish 
between the HCP and FCC lattice structures. Finally, 
the discrepancy between the experimental data and the 
theoretical predictions at neutron energies higher than 
10 meV is accounted by including contributions from in- 
elastic scattering, as explained in the following section. 



FIG. 2: Elastic total cross-sections and experimental data re- 
ported by Seiffert [sJl (solid circles) and our measurements \3^ 
(triangles). Inclusion of the molecular form factor and the 
Debye- Waller factor shifts the polycrystalline cross-section 
prediction from the solid green curve to the solid red curve. 
Elastic cross-sections for both polycrystalline HCP (solid red) 
and FCC (dashed blue) lattice are plotted for comparison. 



the static structure function of single-crystal, which is 
HCP for D2 28|. However, some Raman scattering data 
[23 . [301 indicate an FCC component depending on the 
temperature and the pressure of the solid. Fig. [5] shows 
the calculations together with experimental data. Note 
the molecular form factor and the Debye- Waller factor 
together significantly reduce the elastic scattering ampli- 
tude, as both suppress the neutron scattering amplitude 
associated with large momentum transfers (see Fig. [5] and 
more discussion in Sec lIVp . Note also that the experi- 
mental data do not show the large diffraction peak due 
to scattering off the [Oil] plane. There is evidence that 
the solid D2 grown in a cold cell tends to self-anneal, 
becoming a few single crystals of large size. In these 
experiments, neutrons scatter through diffractions from 
several single crystals of large size oriented at different 
angles. As a result, the total cross-section deviates from 
the prediction in the powder limit described by Eq.(|3]). 

The presented total cross-section is extracted from 
neutron transmission measurements. Excluding the [Oil] 
scattering alone brings the prediction on top of the exper- 
imental data (see the light blue dashed curve in Fig. [2|) . 
In these measurements, the finite area of the neutron de- 
tector would lead to reduced scattering amplitudes, be- 
cause neutrons scattered at low angles are still detected. 
However, the surface coverage of the detector has to be 
unphysically large in order to account for the missing am- 
plitude. This rules out detector geometry as the cause 
of the missing scattering amplitude. On the other hand, 
missing the major [Oil] peak suggests that the crystal is 
grown preferably with the c-axis perpendicular to the di- 
rection of the neutron beam. This would lead to reduced 
Bragg scattering for small momentum transfer centered 



IV. INELASTIC SCATTERING 

The inelastic scattering involving creations and anni- 
hilations of single phonon in Eq. ^ is calculated numer- 
ically for coherent and incoherent processes with corre- 
sponding scattering lengths given by Eqs.(IS]) and PT|) . To 
evaluate coherent scattering in a polycrystalline sample 
beyond the lA prescription, we developed a Monte-Carlo 
algorithm to sample the dispersion curve of each single 
crystallite oriented at random angles, strictly following 
the kinematics of energy and momentum conservation in 
each coherent scattering event [33|- Taking the angular 
average of the momentum transfer vector simplifies the 
orientation-dependent scattering amplitude into a two di- 
mensional map of S(Q, to). The result of one such calcu- 
lation (a high resolution map with 500 x 200 grids in the 
(Q, u) space) is shown in Fig. [31 The calculation includes 
1-phonon and 2-phonon contributions from both the co- 
herent and incoherent process, calculated independently 
for the ortho -^ ortho (J = -^- 0) and the ortho ^• 
para ( J = — > 1) transition. Note that the neutron scat- 
tering resulting in ortho — > para transition involves spin 
flip, and thus is a purely incoherent process. The relative 
weight of the scattering amplitude of the ortho-ortho and 
ortho-para processes is dictated by Eqs.dH]) and (|TT|) . In 
addition, the calculation also includes inelastic scattering 
leading to the rotational excitation Jgi, where incident 
neutrons lose energy in ortho-D2 through inter-molecular 
excitations without activating translational phonons via 
a nuclear recoil. In upscattering, this J = — > 1 process 
is absent in ortho-D2 for neutrons with energy less than 
the Jqi transition energy of 7.1 meV. 

The average algorithm used in polycrystalline samples 
prescribes that, for every given magnitude of momen- 
tum transfer, a random angle is assigned to determine 
the momentum transfer vector. Setting this momentum 
transfer to that of the phonon, the eigenfrequencies and 
eigenvectors are calculated uniquely using the force ten- 
sor matrix. For the HCP structure, the result consists 
of six eigenvalues, each corresponding to a translational 
degree of freedom originated from the two molecules in 
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FIG. 3: A contour plot of the dynamic structure function, 
S(Q,a;), of solid ortlio-D2 at 5 K, plotted for both neutron 
downscattering (positive uj) and neutron upscattering (nega- 
tive uj). To better resolve the details of the phonon modes, 
the upscattering region is plotted on an amplified scale com- 
pared to the downscattering region. The red curves are free 
neutron energy dispersion curves, outlining the allowed kine- 
matics for UCN production and UCN upscattering. The light 
curves form the kinematic bounds for incident neutrons with 
initial energy 30 meV (fine white dots), 10 meV (white solid), 
3 meV (white dashed), and 10 fieV (green dashed-dot). 



a unit cell. The energy transfer bin that contains the 
calculated eigen-frequency is then augmented by a nu- 
merical value of the scattering amplitude estimated using 
Eq.(l2]). For each magnitude of momentum transfer, up to 
10^ different angles are sampled isotropically over 47r. In 
the end, a two-dimensional map of the dynamic structure 
function is evaluated for positive energy transfer, which 
corresponds to neutron downscattering. A second map 
is evaluated independently for negative energy transfer 
(neutron upscattering). The main difference between 
upscattering and downscattering arises from the depen- 
dence on the phonon occupation number n. Upscattering 
occurs when the neutron encounters a phonon, and thus 
the scattering probability scales with the thermal popula- 
tion of phonons. In contrast, downscattering of neutrons 
through phonon creation takes place only when the kine- 
matics of energy and momentum transfer is favorable, as 
determined by the existing modes of phonon excitation. 
The downscattering amplitude scales as (1 -I- n), where 
n weakly enhances the transition probability due to the 
bosonic nature of phonons. 

Next, to evaluate the incoherent contributions, we con- 
struct the density of states by integrating the coherent 
dynamic structure function derived in the previous step, 



1.5 



0.5 



— 0— Yu 








-^^ IVIacFarlane 








Debye Model 




J 


_ 


— Q>Q,rea/2 




Ik 




allQ 




,\ 
















/ R 








/ 1 








'HJ 


\ - 


i 


I ^ 


^w ■ 




/ 


<* 


^^^p 


'1 


I 


• 




t .' 


\^N2ji^i^i0^ 






h 



6 8 10 

Energy (meV) 



12 



FIG. 4: The density of states, Z{ijj), derived from the coherent 
dynamic structure function S'^" (Q, u>) in solid ortho-D2 at 
5 K, with different lower bounds of the integration range in Q. 
The results are compared to the density of states presented in 
Yu [3j] (solid green line) , MacFarlane [3| (solid black line with 
dots) , and that of a simple Debye model (magenta dashed-dot 
line) . 



after correcting all the dynamical factors, over the whole 
accessible range of momentum transfer: 
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(13) 

However, simply integrating over all possible Q values 
leads to a non- vanishing density of states as Q approaches 
zero. This disagrees with the Debye model, which pre- 
scribes the density of states for phonons as cQ^, where c 
is the speed of sound in the solid. To fix this discrepancy, 
we set the lower bound of the integration to half of the 
momentum carried by a free neutron, for the following 
reasons: As illustrated in Fig. [31 the scattering ampli- 
tude peaks beyond the first Brillouin zone. The accessi- 
ble kinematic range of free neutrons in the cold regime 
(with incident energy less than 10 meV) simply cannot 
produce excitations within the first Brillouin zone, except 
for those close to the zone boundaries. In other words, 
the sound speed of phonons exceeds the group velocity of 
the neutron wave-packet, leading to no coupling between 
neutrons and phonons of low momentums. Excluding the 
rather large amplitude for coherent creation of acoustic 
phonons with momentum smaller than half of the free 
neutron momentum (in the first Brillouin zone) restores 
the expected Q dependence at low momentum transfers 
(as shown in Fig|4]). Only the coherent Jqo process is 
included in evaluating Z{uj). 

We have also calculated Z{uj) using the classical ap- 
proach by integrating the derivative of the dispersion 
curve over the hypersurface with constant w, following 
a method described in [S^- The result is consistent 
with the method discussed above. With ^(w), the in- 
coherent scattering law can be calculated using the stan- 



Form Factors in Solid D„ 




FIG. 5: The Debye- Waller factor, e , and the molecular 
form factors for a few lowest rotational transisions. 



dard textbook formula, without the need for the rather 
computation-intensive algorithm developed for the coher- 
ent process. 

For a final consistency check, using the derived density 
of states, we can calculate the Debye- Waller factor using 
the polarization vectors in the HCP structure |25j : 



W{Q) = 3 / dujZ{Lo 



AMuj 



(2nH + l)|g.elL 



(14) 



and examine how well the approximation of the cubic 
symmetry of Eq.Q applies to the HCP structure. We 
found that the Debye- Waller factor of solid D2 at 18 K 
agrees remarkablly well with the simple cubic approxi- 
mation, and that at 5 K is slightly above the value of the 
cubic structure (see FiglS]). For the Q range relevent to 
UCN production (Q < 2.5/A), the agreement is within 
5%. For completeness, the molecular form factors are 
plotted to compare with the Debye- Waller factor. For 
the J = -> transition, the form factor suppresses 
the high momentum transfer, similiar to the effect of the 
Debye- Waller factor. For the J = ^ 1 and J = ^ 2 
transitions, the form factors start at zero and grow with 
increasing Q. Processes involving the J = — > 2 tran- 
sition (E>15 meV) are highly suppressed in magnitudes 
and thus we don't expect it to contribute significaly to 
the cross-section. 

After combining independent results for coherent and 
incoherent cross-section, the major features of the inelas- 
tic scattering are clearly visible in the dynamic structure 
function evaluated with our full model (see Fig. |3]). The 
acoustic phonon branches that extend from elastic Bragg 
peaks into high energy transfers can be clearly identified. 
Only phonons created along the symmetry axes have en- 
ergy dispersion that is periodic in momentum; phonons 
created along other directions have less symmetry, but 
their energy can be calculated analytically using the force 
tensor. We find that the clustering of the scattering am- 
plitude around 5 meV is due to mostly the phonon modes 
in the basal plane. Including scatterings with directions 



out of the basal plane, the phase space increases result- 
ing in large scattering amplitude associated with these 
in-plane phonon modes. The narrow band of large scat- 
tering amplitude around 9 meV is due to the longitudi- 
nal optical branch of phonons associated with scatterings 
along the c-axis. Multi-phonon contribution increases the 
scattering amplitude for energies larger than 10 meV and 

o 

momentum transfer larger than 3/A. 

Overall, the scattering amplitude is significantly sup- 

o 

pressed for momentum transfers larger than 5/A, owing 
to the combined effect of the molecular form factor and 
the Debye- Waller factor. The same suppression due to 
these form factors was observed in the elastic scatter- 
ing as discussed in the previous section. In contrast to 
the interpretations given in |36| , our study shows that 
the large scattering amplitude cannot be simply asso- 
ciated with phonons on the zone-boundaries along the 
symmetry axes, because the phase space for interaction 
with these phonons is very small compared with that 
available for inelastic scattering. Finally, the pure Jqi 
transition (without phonon creation) is visible at a fixed 
energy transfer of 7.1 meV. The molecular form factor 
2ji{Qa/2) dictates the Q dependence of the scattering 
amplitude in the Jqi transition, as it increases from zero 
with increasing Q. By comparison, the dynamic struc- 
ture function calculated in the lA approach (as shown in 
Fig. [1} is grossly over-simplified. 

In order to account for the scattering amplitude be- 
yond the energy of the most energetic phonon available 
(i.e., the transverse optical phonon moving in the [001] 
direction), we have also calculated the 2-phonon contri- 
butions. The coherent dynamic structure function in- 
volving two phonons can be expressed as: 
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where the first term describes the creation of two phonons 
through a single scattering process, the second and the 
third term describe the simultaneous creation of one 
phonon and annihilation of another phonon of different 
energy, and the final term describes the annihilation of 
two phonons in a single scattering event. The combined 
amplitude is smaller than that of the single phonon pro- 
cess. Nevertheless, the two-phonon process enlarges the 
scattering phase space to include those at higher energy 
transfers, and thus neutrons with energies higher than 
the Debye temperature participate in neutron downscat- 
tering through multi-phonon processes. 

In constructing the S{Q,ui) for the two phonon pro- 
cess, we implement the following steps: For a given mag- 
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FIG. 6: The dynamic structure function S'^°''(Q, uj) in solid 
ortho-D2 at 5 K. The crystal orientations in these polycrys- 
talline samples are limited to 10° and 45° relative to the basal 
plane or to the long axis (c-axis). 



nitudc of momentum transfer Q, first throw the dice (i.e., 



the random number generator in the Monte-Carlo code) 
to determine one angle. This random angle together with 
Q determines the vector of momentum transfer. Throw 
the dice again to determine the magnitude of the mo- 
mentum of the first phonon, and again to determine the 
angle of this momentum vector. Once the two vectors 
are specified, the momentum of the second phonon is de- 
termined by the momentum conservation law. Now with 
the momentum vectors specified, we can calculate the en- 
ergy of each of the two phonons and sum them up as the 
total energy transfer. The energy bin corresponding to 
the total energy transfer for this scattering event is incre- 
mented by the amplitude described in Eg. p^ . Note that 
the momentum conservation and energy conservation are 
strictly applied in every scattering event as dictated by 
the delta functions in Eq. p^ . 

To address the effect of preferred directions of crystal 
orientation as indicated by the total cross-section data 
(see Sec lIIip . we have also calculated the scattering laws 
by restricting the scattering angle relative to the inverse 
lattice vectors. A few examples of the resulting dynamic 
structure functions are shown in Fig. |6l To accentuate 
the differences, only the coherent scattering contributions 
are plotted. Comparison with recent measurements |36j 
suggests that the experimental data are best described by 
scatterings confined to within 45° of the basal plane. The 
detailed information contained in these inelastic scatter- 
ing maps strongly suggests that the missing peak in the 
total cross-section (Figl2]) could be explained by a crys- 
tal orientation where the c-axis is aligned preferentially 
perpendicular to direction of neutron propagation. 

With the detailed dynamic structure function, calcula- 
tion of the total cross-section becomes a straightforward 
integration of Eq.([T]) over the appropriate phase space: 
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For neutron downscattering, S'^°"'"{Q, +uj) (upper graph 
in FigIS]) is integrated over the energy transfer w from 
to Eq, where Eq is the energy of incident neutrons. 
For neutron upscattering, S'"^(Q,— w) (lower graph in 
Fig El) is integrated over energy transfer from to cxj. For 
every w, the range of integration on Q is bounded by the 
momentum transfer in forward scattering and backward 
scattering: 
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The range of integration is illustrated by the white lines 
in Fig|3] for several incident neutron energies. Since 
the dynamic structure function derived for polycrys- 
talline target is already angle-averaged, the integration 
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FIG. 7: The calculated total cross-section per D2 molecule in 
solid ortho-D2 at 5 K, plotted in linear scale (upper graph) 
and log scale (lower graph). The results are compared to ex- 
perimental data. The total cross-section is a sum of the elas- 
tic cross-section (green triangles), the downscattering (blue 
x) and the upscattering (red open circles) cross-sections. 



in Eq. (|16p over the azimuthal angle simply yields 2tt. 
In the end; to report the total inelastic scattering cross- 
section per molecule, the result is divided by the number 
of molecules contained in each unit cell, which is two 
in HCP D2 (shown as the pre- factor of 1/2 in Ea. ([TC| ). 
The same normalization is applied to the incoherent com- 
ponent. According to Eq.(|5]), the lattice form factor is 
f^pinc^2 _ 2_ The form factor together with the nor- 
malization brings the incoherent amplitude back to the 
scattering cross-section of one particle. The total cross- 
section is plotted in Fig. [71 For solid ortho-D2 at 5 K, 
the upscattering is negligible for cold neutrons over the 
plotted energy range from 0.1 to 20 meV, whereas the 
downscattering amplitude becomes appreciable for inci- 
dent neutron energy larger than 5 meV. Excluding the 
missing [Oil] Bragg peak ranging from 3 to 6 meV, the 
inclusion of contributions from inelastic scattering pro- 
duces reasonable agreements to the total cross-section 
data derived from simple transmission measurements. 



V. UCN CROSS-SECTIONS 

The cross-sections for UCN production and UCN up- 
scattering can be evaluated simply by finding the value 
of S{Q,uj = Efree) whcrc Efree — h^Q^/2m, which is 
the dispersion curve of free neutron (the parabolic curve 
shown in Fig. [S]). Here, the integration of Eq. lfTH)) sim- 
plifies because the energy scale of UCN (< 350 neV) 
is several orders of magnitude smaller than the energy 
scale of cold neutrons (^1 meV). For inelastic scatter- 
ing events leading to UCN productions, the momentum 
transfer Q ~ ki — kucn equals the initial momentum of 
the neutron ki due to the negligible fc„cn- Similarly, the 
energy transfer lo = Ei — Eucn ~ Ei. In estimating the 
UCN production cross-section, the integration of the dou- 
ble differential cross-section Eq.(IT|) becomes: 
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A change of integration variable dEucnd^f — 
dEucnd(j)dii = d(j)Y^ — dujdQ, where /x = cos{6), further 
simplifies the above integral to 
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For scattering with the final state of interest, i.e., UCN 
with energy Eucn, the range of integration is a narrow 
band around the parabola of the free neutron, Q* = 

ki, with w* 
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Ej, and thus the above in- 



tegral can be reduced to: 
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(18) 



given that S(Q, ui) is a smooth function. Next, to esti- 
mate the production of UCN that can be confined in the 
experimental apparatus, in which the Fermi-potential of 
the container material sets the maximum energy of UCN 
(Vf = Eucn), an additional step is required to integrate 
the UCN energy from to Vp = ffkl^^/2m: 
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As a result, the cross-section of UCN production for 
neutrons with incident energy Ei is 
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(19) 
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UCN Production in Solid ortho-D (Detailed Model) 
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FIG. 8: Upper: The UCN Production cross-section bro- 
ken down into full different processes included in the full 
model, including coherent and incoherent excitations of 1 and 
2 phonons, coupled to the two lowest rotational state tran- 
sitions, i.e., Joo and Joi. The Joi process without phonon 
excitations also contributes to the UCN production. Lower: 
The prediction using the full model is compared to that based 
on the lA and the simple Debye model, along with the ex- 
perimental data. Energy spectra of cold neutron flux with 
Maxwell-Boltzmann distribution at 40 K and 70 K are plotted 
to illustrated the degree of overlap with the energy-dependent 
cross-sections. 



to produced storable UCN with Eiji^^n from to E^cn = 
Vf. 

Two sets of experimental data on UCN production are 
compared to the predictions of the full model. First, note 
that the production cross-section reported in [9| is nor- 
malized to each atom and their model does not include 
either the spin statistics or the molecular form factor of 
molecular D2; the cross-section independently reported 
hi [36,..37|| should be corrected by a factor of two to prop- 
erly account for the range of integration on the allowed 
momentum transfer through — fc„cn to +kucn- To allow 
for direct comparisons, we multiply both data sets by 
a factor of two without any further corrections. Over- 
all, the full model gives fair agreement with both the 
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FIG. 9: The optimization of UCN Production by varying the 
incident cold neutron spectrum. The optimized cold neutron 
flux is a 33 K Maxwellian for I A treatment. The optimized 
cold neutron spectrum is a 40 K Maxwellian evaluated with 
the full model. 



results of indirect UCN production extracted from the 
data of cold neutron scattering [33| and the results of di- 
rect UCN production [3(as shown in Fig. [8]). Around the 
peak production at 6 meV, the full model gives a strik- 
ing agreement with the experimental data, whereas the 
lA approach using a somewhat realistic density of states 
misses the peak by 1 meV. The simple Debye model fails 
to capture any details of energy dependence. The Jqi 
transition at 7.1 meV is a delta function in the model, 
but the finite energy resolution of the scattering instru- 
ment integrating around the transition energy reduces 
the amplitude. In spite of it, the data suggests that the 
transition energy is slightly higher than 7.1 meV, as used 
in the calculation. The second peak around 9 meV pre- 
dicted by the full model is absent in the Miinich data [331 ^ 
however, is present in the PSI data. The disappearance 
of the 9 meV peak could be a result of the preferred crys- 
tal orientation as discussed in Sec lIVI By restricting the 
scattering angle around the basal plane, the scattering 
intensity of the 9 mcV peak can be adjusted to a smaller 
value (as shown in Fig. |6]). For low scattering intensities 
at energy range below 4 meV and higher than 11 meV, 
the full model predicts production cross-sections smaller 
than the reported experimental data. This might simply 
due to the artifact of insufficient background subtraction 
in the experimental data. In particular, the subtraction 
of non-vanishing tails of the large elastic peaks is diffi- 
cult to carry out. As for the high energy end, multiple 
scatterings couple to neutrons of energies higher than 
the Debye temperature, and thus enhance the scatter- 
ing amplitudes. Neither of the data points is corrected 
for effects of multiple scattering. The very different di- 
mension of the target used in these two experiments also 
hints that the origin of the high energy excess might be 
multiple scattering in origin. 

UCN production can be optimized by varying the en- 
ergy spectrum of the incident cold neutrons. Integrat- 
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Upscattering of UCN (500A) in solid o-D 
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FIG. 10: The evolution of the differential cross-section 
Cpa{Eucn -^ Ef)/dEfdQ for UCN upscattering with increas- 
ing source temperature. Results of the full model (solid 
curves) and the lA predictions (magenta dashed curves) are 
plotted as a function of the final energy of upscattered neu- 
trons. The inset plots the energy-integrated contributions of 
the coherent and incoherent process in the full model, as a 
function of the source temperature. 



ing the product of the energy-dependent cross-section for 
UCN production and the incident neutron flux over the 
energy spectrum yields a UCN production cross-section, 
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that depends on the temperature of incident cold neu- 
trons. Here we assume that the cold neutron flux has a 
thcrmalized energy spectrum with temperature Tc„. As 
shown in Fig.[9l the UCN production is the greatest when 
coupled to a flux with a 40 K Maxwell-Boltzmann energy 
spectrum, whereas the optimized cold neutron spectrum 
predicted by the lA is 33 K with a cross-section 32% 
larger than that predicted by the full model. This dif- 
ference comes from the fact that the lA approach over- 
estimates the contributions below 5 meV, where the free 
neutron parabola does not intersect with most phonon 
branches at low Q (see Sec. lIVp . Even though the differ- 
ential cross-section for UCN production peaks at 6 meV, 
the overlap with a flux of cold neutrons with a 70 K spec- 
trum does not result in more UCN production, because 
many neutrons in the Boltzmann distribution are spread 
out over the higher energy range resulting in a reduction 
of peaked flux. Using a colder flux of neutrons, a higher 
percentage of neutrons are directly under the production 
peak. 

The achievable UCN density does not only depend on 
the production cross-section, but also scales with the life- 
time of UCN, which sets the limit on the time duration 
during which the UCN density can accumulate without 
loss inside the source. Upscattering is a source of loss that 



needs to be controlled. Following the same approach for 
UCN production, the cross-section for UCN upscattering 
is evaluated by finding the S{Q, — w) along the free neu- 
tron parabola as illustrated in the lower graph in Fig. [31 
Note that the upscattering amplitude is directly propor- 
tional to the occupation number of phonons, which obey 
the Bose-Einstein statistics. The occupation number 
peaks at low w and falls off in magnitude as w increases. 
The amplitude of UCN upscattering depends even more 
strongly on the physical presence of these phonons than 
that of the downscattcring. For small momentum trans- 
fers within the first Brillouin zone, the free neutrons do 
not intersect with any branches of acoustic phonon (see 
the lower graph in Figl3|). Without direct coupling to 
the acoustic phonons, there is no coherent process that 
upscattcrs UCN though phonon annihilation. For solids 
at low temperatures, the non-zero upscattering ampli- 
tude comes from the remaining incoherent 1-phonon and 
multi-phonon processes (see Fig. [T0|) . As the tempera- 
ture increases, the contribution from coherent scattering 
increases with increasing thermal population of phonons 
beyond the first Brillouin zone (shown in the inset plot 
in Fig. I10|) . and the coherent phonon annihilation take 
place. Fig. [10] presents the evolution of the differential 
cross-section of UCN upscattering with increasing tem- 
perature from 5 K to 18 K. 

The total cross-section for UCN upscattering is calcu- 
lated by integrating S'^p{Q,—uj) (lower graph in Fig[3]) 
following: 
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Upon upscattering, UCN scatters into a narrow band of 
phase space with the final energy and momentum around 
the free neutron parabola (as illustrated in the lower 
graph in Fig. [3|). The integration simplifies into: 
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In Fig. [TT] plots upscattering cross-section for UCN with 

o 

a wavelength of 500 A. For UCN with a different wave- 
length, the upscattering is rescaled linearly with its initial 
wavelength. The full model predicts a UCN upscatter- 
ing cross-section 2^4 times smaller than that estimated 
by lA (sec Fig. fTTj) depending on the temperature. The 
excess in the lA prediction is due to the use of the den- 
sity of states that over-estimates the number of modes in 
small energy transfers within the first Brillouin zone (see 
details shown in Fig.fTOj). This correction is by no means 
small. With a reduced upscattering cross-section, UCN 
can live longer inside the source at higher temperatures. 
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UCN (500 A) upscattering cross-section in solid ortho-D 
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FIG. 11: UCN upscattering cross-section (per D2 molecule) 
in solid ortho-D2 as a function of the source temperature. 
Both the full model (blue line with x) and the lA predic- 
tions (red line with open circles) are plotted. The nuclear 
absorption (green line with solid circles) is independent of 
the source temperature. The intersection points between the 
upscattering cross-section and the nuclear absorption cross- 
section (square box) sets the saturation temperature of UCN 
yield. The ratio between the lA and full model prediction is 
plotted (dashed line) and presented along the y-axis on the 
right. 



The ultimate limit on the UCN lifetime in solid D2 
comes from the small, yet non- vanishing neutron absorp- 
tion of each deuteron nuclide. In a source where the es- 
cape time of UCN is comparable to (or larger than) the 
absorption time, the UCN yield saturates when the up- 
scattering loss is reduced to the same level as the nuclear 
absorption loss. Reduction of the upscattering by fur- 
ther cooling would not increase the UCN yield of UCN 
by more than a factor of two. In practice, the gain is 
smaller than two because of the additional sources of 
loss, such as the hyd rog en contamination and the para- 
D2 contamination [8J, |39[ . The condition for saturation in 
UCN output sets the practical operational temperature 
of the source, and determines the requirements on cryo- 
genics engineering. The direct consequence of reduced 
upscattering cross-section in solid ortho-D2 is that the 
predicted saturation temperature is ^ 7.5 K, which is 2 
degrees higher than that predicted previously using the 
lA. The higher operating temperature makes the imple- 
mentation of cryogenic source somewhat less challenging 
as most refrigerator systems have higher cooling power 
at elevated temperatures. 

For a practical sample of converted D2 gas containing 
a residual 3% of para-D2, the total loss cross-section is 
increased by 1 b, because the temperature-independent 
upscattering cross-section of pure para-D2 is 31 b Q. In 
this case, the full model predicts the saturation temper- 
ature for UCN production to be 11 K, and the I A model 
predicts it to be 8 K. Using un-convertcd normal D2 gas 
with 33% para-D2, the loss cross-section is about 10 b. 
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UCN (500 A) upscattering cross-section in soiid ortlio-D 
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FIG. 12: Upscattering cross-sections (for 500 A UCN) mea- 
sured in a UCN experiment [32|, compared to predictions from 
the full model and the IA. The vertical line indicates the triple 
point of solid D2 at temperature of 18.678 K, separating the 
solid and liquid phases. 



According to the upscattering cross-sections presented in 
Fig. [11] the UCN yield would reach saturation at temper- 
atures as high as 18 K, and no superthermal gain would 
be observed by cooling the UCN converter. 

Finally, the updated upscattering cross-section pre- 
dicted by the full model of S{Q,uj) for 0-D2 is applied 
to understand the UCN production data from our recent 
measurements using solid D2 |32|]. The Monte-Carlo sim- 
ulation using the upscattering cross-sections calculated 
using the IA approach fails to reproduce the tempera- 
ture dependence measured in the experiment. Details of 
the Monte-Carlo simulation can be found in JS^]. We 
also investigated other effects that lead to increased elas- 
tic scattering inside the source, however, we find that the 
higher saturation temperature can only be explained by 
reduced loss. The coherent scattering leading to smaller 
upscattering cross-sections predicted using our full model 
is a very likely candidate to provide the required modifi- 
cation. 

We then extract upscattering cross-sections using the 
same set of experimental data reported in [32| . The up- 
scattering cross-sections of solid D2 are derived from an- 
other Monte-Carlo study, in which the upscattering cross- 
section is varied until the simulated result on UCN yield 
best fits the experimental data for each different temper- 
ature. Results of this study arc shown in Fig. [121 com- 
pared to theoretical predictions based on the full model 
and the IA. The Monte-Carlo code includes the full range 
of UCN energy spectrum from up to 1 ^eV. To facil- 
itate direct comparison with previous calculations, the 
extracted upscattering cross-section are for UCN with a 
wavelength of 500 A. Below the triple point, the upscat- 
tring cross-sections agree quite well with the full model 
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FIG. 13: The temperature dependence of our FP12 exper- 
iment data reported in [34I is compared to similar experi- 
ments (Mainzgil, LANSCE UCNA[3i|, PNPl[43|, PSljiJl). 
All data are scaled to unity at 5 K. GEANT4 simulations 
using upscattering cross-sections calculated by the full model 
(solid curve) is compared to that using the I A (dashed curve). 



prediction. Above the triple point, the mechanism for 
UCN upscattering (in the liquid phase) is significantly 
different from that in the solid, and none of the models 
discussed in this paper apply. Data points of small up- 
scattering cross-sections fluctuates at low temperatures 
due to low statistics of UCN signal in general. At tem- 
peratures higher than 10 K, where the UCN upscattering 
cross-section is larger than 1 b, the incoherent model is 
excluded at the 2cr level and higher. 

We also surveyed experimental results on UCN pro- 
duction using solid deuterium reported independently by 
various research groups. Data sets are plotted in Fig. [13] 
together with the results of simulations. While the re- 
sults of simulation using the updated upscattering cross- 
section agree quite well with our data and the data set 
measured at PSI, the data sets from LANL and Mainz 
groups show a much steeper temperature dependence for 
temperatures higher than 10 K. The difference comes 
from the different source configurations. In these two 
experiments, solid D2 was condensed from vapor at the 
end of the UCN guide, which was cooled below the so- 



lidification temperature of D2 . The source was designed 
to reduce the transmission loss by eliminating the vac- 
uum window that is typically installed to contain the 
volatile D2. The windowless source worked quite well at 
low temperatures, however, for temperatures higher than 
10 K, the entire UCN guide was filled with D2 gas at the 
saturated vapor pressure. The large upscattering cross- 
section from the D2 vapor, together with a large volume 
where the vapor permeated, resulted in a steeper temper- 
ature dependence than the simple prediction where UCN 
upscatters through absorption of phonon in the confined 
region of the source. Both PSI and our experiment used 
D2 contained in a target cell, in which the vapor region 
is considerably smaller. None of the data sets show a 
low saturation temperature, as predicted by the lA. In- 
stead, the temperature dependence of experimental data 
points agrees better with results of the simulation using 
upscattering cross-sections predicted by the full model. 



VI. CONCLUSION 

In this paper, we have presented calculations of the 
neutron cross-sections in polycrystalline ortho-D2. All 
known physics on spin statistics, rotational excitations, 
elastic Bragg scattering, coherent inelastic, and incoher- 
ent inelastic scattering have been implemented in our 
full model for the first time. The presented work ad- 
dresses the shortcomings of the widely used lA in esti- 
mating neutron cross-sections at energies at thermal and 
cold regime. For applications on modern neutron sources 
for UCN production using solid D2 and other materials 
with large coherent scattering lengths, such as helium, 
oxygen and nitrogen, the interference between scattered 
neutron wavefunctions across all lattice sites is important 
and should not be ignored. For the case of solid ortho-D2, 
we have shown that UCN cross-sections are significantly 
altered under the new treatment. In particular, the UCN 
upscattering cross-section is a factor of 2^4 smaller than 
previous predictions throughout the temperature range 
of interest to UCN production. The same modifications 
will impact designs of sources (beyond UCN sources) that 
aim to produce long wavelength neutrons using materials 
other than hydrogen. The Monte-Carlo algorithm devel- 
oped in this work can be applied to all materials with sig- 
nificant coherent scattering components, provided that 
the energy dispersion of the quasi-particles is known. 
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